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Baryon fractions in clusters of galaxies: evidence against a 
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ABSTRACT 

The MillenniuiTL Gas project aims to undertake smoothed-particle hydrodynamic 
resimulations of the Mihenniuni Simulation, providing many hundred massive galaxy 
clusters for comparison with X-ray surveys (170 clusters with kT^i > 3keV). This 
paper looks at the hot gas and stellar fractions of clusters in simulations with differ- 
ent physical heating mechanisms. These fail to reproduce cool-core systems but are 
successful in matching the hot gas profiles of non-cool-core clusters. Although there is 
immense scatter in the observational data, the simulated clusters broadly match the 
integrated gas fractions within r^oo- In line with previous work, however, they fare 
much less well when compared to the stellar fractions, having a dependence on cluster 
mass that is much weaker than is observed. The evolution with redshift of the hot 
gas fraction is much larger in the simulation with early preheating than in one with 
continual feedback; observations favour the latter model. The strong dependence of 
hot gas fraction on cluster physics limits its use as a probe of cosmological parameters. 
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1 INTRODUCTION 

The hot gas fraction of clusters of galaxies was first used as a 
co smological probe b y lAUen et al . (2002), and later refined 
in I Allen et ail (|2004l ') and lAllen et all (|200i ). These papers 
show quite conclusively that the gas fraction can be used to 
derive cosmological parameters that are in agreement with 
the concordance ACDM cosmology. A similar result was ob- 
tained by LaRoquc ct al. (2006) in a joint X-ray, Sunyaev- 
Zel'dovich analysis. Subsequently, a more sophisticated anal- 
ysis that properly takes into account selection effects, and 
combines mass. X-ray luminosity and temperature observa- 
tions of 238 c lusters at z Sj 0.5, reached similar conc l usions 
iJMantz et a l.l l2009 a.b). as did a study bylEttorTeFaLl (|2009l ) 

of 60 clusters extending to 2: '^ 1.3. 

However, as has been pointed out by ISadat et al.l 

(|2005l ). the above conclusion relies very heavily upon the 
assumption that the gas fraction in the clusters used in the 
study is independent of both mass and redshift. We show in 
this paper that using different models for entropy generation 
in the intracluster medium (ICM) can lead to variations in 
hot gas fractions in simulated clusters that are at least as 
great as those that one obtains by using an incorrect cosmol- 
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ogy in the observational data analysis. We argue, therefore, 
that, at present, it is more useful to fix the cosmology to 
the concordance value and to use the data to constrain clus- 
ter physics. By doing so we conclude that the data favour 
a model of continual energy injection into the ICM from 
galaxies rather than a widespread preheating episode at high 
redshift. 

In the future, once the physical models of the ICM be- 
come more refined, and with the large statistical samples of 
clusters generated by eROSITA, it should be possible to do 
a combined analysis that constrains both the cluster physics 
and the cosmology simultaneously. 

In Section [2] we describe the numerical method that 
we use, the different feedback schemes, and our method of 
cluster identification. Section [3] describes our results first on 
gas fraction profiles, then scaling relations, and finally the 
evolution of each of these. The conclusions of the paper are 
summarised in Section [4] 



2 METHOD 

Here we present an overview of ou r numerical scheme. The 
method is described at length in IShort fc Thomad (|2009l ) 
and its particula r application to clusters of galaxies in 
IShort et all (|2010l ). 
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2.1 Simulations 

We present results from three simulations taken from the 
Millennium Gas Project, the basic objective of which is 
to add ^as to the da rk matter-only Millennium Simulation 
l|Springel et al.1 12003 ) . Each simulation incorporates a dif- 
ferent model of the baryonic physics, so that we can assess 
the impact of varying physical assumptions on the thermal 
history of the ICM. 

In the first Millennium Gas run, the intracluster gas 
is heated solely by gravitational processes. We refer to this 
simulation as the GO (Gravitation Only) run. Although this 
run does not include gas cooling or heating from astrophys- 
ical sources such as supernovae (SNe) and active galactic 
nuclei (AGN), it is useful as a base model, enabling us to 
determine exactly which cluster properties are affected by 
astrophysical processes beyond gravitational heating. Given 
that the only source of gas entropy changes in the GO 
run is gravity, then we would expect a self-similar clus- 
ter population to be formed. This is gene rally found to 
be the case in non- radiative simulations f e.g. iNavarro et al.l 
19951; lEkeet al.|[l99 8: Voit ot al. 2005; Asc asibar et al.ll2006l : 



Muanwong et al.ll2bog : IStanek et al.ll20ia ). 



The second Millennium Gas simulation also includes 
high-redshift pr e heati ng (first posited by iKaisen Il99ll : 
lEvrard fc Henrvl 1 199 il l and radiative cooling. We name 
this simulation the PC (Preheating plus Cooling) run. 
Preheating raises the entropy of the ICM before gravita- 
tional collapse, preventing gas from reaching high densi- 
ties in central cluster regions and thus reducing its X- 
ray emissivity. This effect is greater in lower-mass sys- 
tems, breaking the self-similarity of the clust er scaling re- 
lations in a way that resembles o b servations |Bialek et al.l 
200ll; iBrighenti fc Mathews! l200l]: iMuanwong et all l200a 



Borgani et all |2002| : iTornatore et all l2003l : iBorgani et al] 



200, 



The simple model of pre heating employed in the PC 
simulation is similar to that of lBorgani et al.l (|2002l ) . Briefly, 
the entropy of every particle is raised to 200 keV cm^ at 
z = 4, thus creating an entropy 'floor' (note that a par- 
ticle's entropy is not changed if it already has a value in 
excess of this at z = 4). In addition to preheating, there 
is also radiative cool i ng ba sed on the cooling function of 
[Sutherland fc Dopital (ll_993'), assuming a fixed metallicity of 
0.3 Zq (a good approximation to the mean me tallicity of 
the ICM out to at least z = l: lTozzi et ai]|2003l ). Once the 
temperature of a gas particle drops below 2 x 10* K, the 
hydrogen density exceeds pn = 4.2 x lO^'^^gcm^'^ and the 
density contrast is greater than 100, then it is converted to a 
coUisionless star particle. However, the preheating is so ex- 
treme that star formation is effectively terminated at z — 4, 
so that less than 2% of the baryonic matter is locked-up in 
stars at 2 = 0. 

Although there is considerable evidence that non- 
gravitational heating of the ICM indeed occurs mainly 
at hi gh redshift (e.g. lEisenhardt et al.l l2008l : IWeiner et al.l 
I2OO9I ). the preheating scenario is clearly a gross simplifica- 
tion of the complex interplay between star formation, black 
hole growth and associated feedback. Despite this, preheat- 
ing does provide a useful effective model for the effects of 
non-gravitational heating. In particular, the Millennium Gas 
PC run can reproduce several key observational properties of 



the lo w-redshift cluster p opulation, including halo gas frac- 
tions JStanek et al.ll2010l ). However, the model fails to ac- 
count for the observed scatter about the mean relations, par- 
ticularly on group scales, and generates over-large isentropic 
core s in low-mass systems as compared to ob servational data 
(e.g. iPonman et alll2003l : iPratt et al.ll2006l '). 

Finally, we use a recent addition to the Millennium Gas 
suite in which feedback is directly tied to galaxy formation, 
rather than assuming some ad hoc injection of energy at 
high-redshift. We term this the FO (Feedback Only) run to 
emphasise the fact that the model currently does not include 
ra diative cooling. The mo del we use is the hybrid scheme 
of IShort fc Thomad (|2009l ). where a semi-analytic model is 
used to calculate the energy transferred to the intracluster 
gas by SNe and AGN. An immediate benefit of this approach 
is that feedback is guaranteed to originate from a galaxy 
population whose observational properties agree well with 
those of real galax;ies. This is generally not the case in fully 
self-consistent hydrodynamical simulations that include ra- 
diative cooling and stellar feedback because too much gas 
cools out of th e hot pha se, leading to excessive star forma- 
tion (e.g. Borg ani et al.l l2004: Kay et al. 2007aj). It is widely 
thought that additional heating from AGN is the natural so- 
lution to t his overcooling probl em. Indeed, iMcCarthv et al.l 
(|2009|) and I Fabian et al.l (|2010l ) have demonstrated that in- 
cluding AGN feedback in hydrodynamical simulations can 
successfully balance radiative cooling in galaxy groups. How- 
ever, the stellar fraction is still found to be 2 to 3 times larger 
than observed in massive clusters. 

For reasons of computational efficiency, the FO run was 
not undertaken over the entire Millennium Simulation vol- 
ume. Instead, we resimulated a sample of several hundred 
galajcy groups and clusters. Each run came in three dis- 
tinct stages: a dark matter-only resimulation of each re- 
gion containing a cluster from our sample; semi-analytic 
galaxy catalogues built on the halo merger trees of these res- 
imulations; and hydrodynamical resimulations of the same 
regions to track the energy injection from model galax- 
ies. The m erger trees were built using the procedure of 
ISpringel et al. (2005j) and we used the Munich L-Galaxies 
se mi-analytic model wi t h the same parameters as described 
in iDe Lucia fc Blaizoj ((20071). Energy is injected into the 
ICM from both SNe an d AGN following the prescription of 
Short fc Thomas (^2009|). The procedure is described in full 
in Short et al. (201Q ) 

By coupling a SAM to a hydrodynamical simulation, 
IShort fc Thomaj (|2009l ') showed that their hybrid feedback 
model could reproduce the observed mean Lx-Tx relation 
for groups and poor clusters at 2 = 0, but only if there 
was a large energy input into the ICM from AGN over the 
entire formation history of halos. The AGN heating is effi- 
cient at driving X-ray emitting gas from the central regions 
of low-mass halos, reducing their luminosity and steepening 
the Lx-Tx relation as desired. Unlike the simple preheating 
scenario, their model was also able to account for some of the 
scatter about the mean relation seen for temperatures T < 3 
keV, attributable to the varied merger histories of groups. 
In addition, the gas fractions of their simulated groups and 
poor clusters were found to broadly agree with observational 
data, rapidly declining at low temperatures and exhibiting 
a comparable amount of scatter. 

The main limitation of both the PC and the FO simu- 
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lations is that neither can reproduce the low entropy found 
in the centres of cool-core clusters. For the FO run this is 
because cooling is not incorporated in the hydrodynamical 
simulations; whereas in the PC run preheating expels gas 
from cluster cores at high redshift, limiting the subsequent 
cooling. 

The cosmological model adopted in all three Millennium 
Gas simulations is a spatially-flat ACDM model with pa- 
rameters Q.^fi = 0.25, fib.o = 0.045, r^A.o = 0.75, h = 0.73, 
Us — 1 and (Ts.o = 0.9. Here f^m.o, fib.o and r2A,o are the 
total matter, baryon and dark energy density parameters, 
respectively, h is the Hubble parameter in units of 100 km 
s"'^ Mpc""*^, Us is the spectral index of primordial density 
perturbations, and crg.o is the rms linear density fluctuation 
within a sphere of radius Sh~^ Mpc. The subscript signifies 
the value of a quantity at the present day. These cosmolog- 
ical parameters are the same as those used in the original 
Millennium simulation and are consistent with a combined 
analysis of the first-ye ar Wilkinson Micro wave Anisotropy 
Probe (WMAP) data (|Spergel et all 120031 ') ar id data from 
the T wo- degree- Field Galaxy Redshift Survey ()Colless et al.l 
I2OOII ). However, there is some tension between the cho- 
sen parameter values, particularly rig and qg.p, and those 
deriv ed from the seven-year WMAP data (JKomatsu et al.l 
I2OI0I ). More significantly for this paper, the mean value of 
the baryon density, /b = 0.18, is higher than the WMAP 
7- year value of /b = 0.167. 



2.2 Cluster catalogues 

Cluster catalogues are generated at several redshifts for the 
three Millennium G as simulations using a pr ocedure similar 
to that employed bv lMuanwong et al.l ([2002). Essentially, a 
friend-of-friends algorthm is used to identify peaks in the 
density field and then spheres are grown around these peaks 
until they enclose regions of a given overdensity. 

We define overdensity. A, with respect to the critical 
density, pc, at any given redshift. 



P(< 



2GM{< r) 

r^m ' 



(1) 



where p is the mean density within radius r, M is the mass 
within this region and G is the gravitational constant. We 
express our masses in units of /i~^Mq, rescaling observa- 
tional data to this system when required. 

Low-mass clusters are more affected by non- 
gravitational heating processes than high-mass ones. 
Observationally, temperature is often used as a proxy for 
mass as the two are strongly correlated. For that reason, 
we often divide the clusters into bins according to their 
spectroscopic-like temperature, T^i, defined as 



Tsi = 



JVr3dv_ 
/p2r-tdF' 



(2) 



where p is the gas density, T the physical tempe r ature , 
and the integral runs over volume. Maz zotta et al.l (I2OOJ) 
have shown Tbi to be a good approximation to the temper- 
ature recovered by X-ray spectral analysis software in the 
bremsstrahlung regime. 

The scaling relations for the GO and PC runs can be 
dominated by small objects. For this reason, we remove 



many clusters from our sample such that the remaining clus- 
ters are distributed evenly in log(M2oo), with a lower mass- 
hmit M2500 > 1-73 x 10^^ h~^ M© that corresponds to 1000 
particles each of gas and dark matter within -R2500- For the 
FO run we resimulate all clusters with kT^i > 3keV and 
a selection of clusters below this temperature, again cho- 
sen evenly in log(M2oo). This has a much higher mass res- 
olution and so the lower mass limit in this case, M200 > 
1.2 X 10^"^ h~^ Mq, was instead fixed by the total number 
of clusters that we wish to simulate. When plotting as a 
function of mass at an overdensity other than the one used 
in the selection procedure, there will not be a clean lower 
mass-limit; however this is a very minor effect that does not 
lead to any significant bias in our results. 



3 RESULTS 

In this section we first discuss the radial profiles of the hot 
gas fractions of clusters at the present day. We then char- 
acterise the dependence of the hot gas fraction upon cluster 
mass, and investigate the scatter about that mean relation. 
Finally, we look at the variation of the hot gas fraction with 
redshift. 



3.1 Profiles 

3.1.1 Differential hot gas profiles 

Figure [T] shows the differential gas mass fraction profiles for 
clusters in three temperature ranges, kT^i > 5 keV (upper, 
green regions), 2.5 keV < fcTsi < 5 keV (middle, cyan regions) 
and fcTsi <2.5keV (lower, magneta regions). The curves are 
plotted out to beyond rsoo (the radius at which A = 500), 
which is the region accessible to X-ray observations. 

For the GO run the gas fraction plateaus at a value of 
0.16-0.17 at about 0.3 7-500, although there is a very slow 
increase at larger radii (beyond the right-hand edge of the 
plot). This is less than the global baryon fraction of 0.18: 
conversion of kinetic energy into heat, together with contin- 
ual stirring of the gas by the motion of dark matter struc- 
tures, allows the g as to pick up e nergy at the expense of the 
dark matter (e.g. IPearce et al.lll994i) . This is particularly 
evident in the cluster cores: for the largest clusters, 7-500 is 
of order 1 /i~^Mpc; therefore the drop in baryon fraction in 
the core of the clusters occurs on a scale significantly larger 
than the force softening (25/i~^kpc). 

In the PC and FO runs, gas has been expelled from the 
cluster cores and pushed to larger radii. The effect is more 
pronounced at lower masses: thus the gas profiles of clusters 
with fcTsi < 2.5 keV are still steeply rising at 7-500, while 
for fcTai > 5keV the profiles are approximately constant 
beyond this radius. The inconstancy of the gas fraction is 
both a nuisance, requiring careful calibration before we can 
use clusters as cosmological probes, and a useful test of any 
model of entropy generation in the ICM. 



3.1.2 Comparison with observations 

Rather than plotting differential profiles, observational stud- 
ies tend to report the cumulative gas fractions, averaged 
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Figure 1. DiflFerential hot gas fraction profiles. The sohd lines 
show the mean relations and the shaded region the l-sigma 
scatter. The upper, green regions correspond to clusters with 
spectroscopic-like temperatures above 5 keV; the middle, cyan re- 
gions to 2.5— 5keV and the lower, magenta regions to clusters in 
the range 1— 2.5keV. Profiles are only plotted for radii greater 
than the gravitational softening length. The dashed line in each 
case shows the cosmic mean. 



within some radius. Figure [2] shows the cumulative gas frac- 
tion profiles in the FQ run compared to several different ob- 
servational studie s (|Vikhlinin et al.ll2006l : lAllen et al.ll2008l : 
iPratt et al.ll2010l ). We have only plotted simulated clusters 
in temperature ranges corresponding to those of the various 
observational data sets. To save space, we do not show the 
results from the PC run — these are similar. 

All three observational studies plot the gas fractions us- 
ing a d ifferent ordinate. The top panel, from lVikhlinin et al.l 
l|2006r) . uses overdensity relative to the critical density. They 
looked at 13 Chandra clusters with a range of temperatures 
upwards of about 1.5 keV for which the data extend out to 
large radii. The simulated and observed clusters agree at the 
outer limit of the data, but the former fall more rapidly as 
one moves into the cluster centre. 
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Figure 2. Comparison of the cu mulative hot gas fractio n pro- 
files in the FO run with those of IVikhlinin et al.l ||2006|. upper 
panel), lAUen et all ||2008| . middle panel) and iPratt et all 1I2OIOI . 
lower panel). In each case the coloured bands refer to the 1- 
sigma spread of profiles seen in the simulations: upper, green 
(fcTai > 5keV); middle, cyan (2.5keV< fcTgi <5keV); lower, ma- 
genta (lkeV< fcTj,! < 2.5 keV). The profiles of simulated clusters 
are only plotted for radii greater than the softening length. The 
observational data are shown by black lines with, in the lower 
panel, solid and dashed lines corresponding to non-cool-core and 
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The middle panel shows data from lAUen et alj l|200a ). 
They again study Chandra clusters, but they focus on the in- 
ner regions of 42 hot (fcTsi > 5 keV) systems. As can clearly 
be seen, the simulated clusters lie well below the observa- 
tions within r25oo- 

Both the above studies focus on bright, relax;ed systems 
(those that are likely to be labelled co ol core, or CC). By way 
of contrast, the REXCESS survey JBohringer et all I2OO7I : 
IPratt et alJl201Cl ). shown in the lower panel of Figure [51 is 
sample of 33 nearby galaxy clusters from XMM-Newton, se- 
lected so as to sample a broad range of luminosities and 
with no bias towards any morphological type. The temper- 
ature range here is 2-9 keV with the more massive clusters 
lying towards the upper edge of the observed band, and the 
least massive ones towards the bottom. Here the simula- 
tions are much more successful in reproducing the observed 
profiles, providing a fair match to the non-cool-core (NCC) 
population out to the limit of the observations. They fail to 
reproduce CC clusters (those with fiattened baryon fraction 
profiles in their inner regions), however these are much less 
frequent than in the relaxed samples. 

We conclude that our simulated clusters provide a fair 
match to the hot gas fractions in typical NCC clusters, but 
fail to reproduce the higher gas fractions seen in the central 
regions of the brighter, CC clusters. 

The results for our adiabatic halos agree with 
other previous simulati o ns undertaken with SPH 
(e.g. iKravtsov et~aLl l2005l : lEttori et al.l I2OO6I : ICrain et all 
[20071) . Direct comparison of the other runs is more difficult 
because we use different feedback models. 

Qualitatively, we see a similar behaviour in the profiles 
of the hot gas fraction to previous work, but the baryon 
fraction profile is very different in our EG run because of 
the much reduced stellar fraction. The integrated baryon 
fractions within rsoo are considered in Section [3.2.31 



3.2 Scaling relations 

3.2.1 Cumulative hot gas fractions 

Figure [31 shows the cumulative gas fraction within a radius 
of rsoo as a function of total mass. The GO points are con- 
sistent with a constant value of 0.162, slightly smaller than 
the universal mean of 0.18. By contrast, both the PC and 
the EG runs have hot gas fractions that are strong functions 
of mass, because the feedback processes are more effective 
in lower-mass clusters and evacuate more of the gas. 

Figure [31 contrasts the hot gas fractions for the EG run 
within three different radii corresponding to enclosed over- 
densities of A = 200, 500, and 2500 (the PC run gives 
similar results). In each case the mass has been measured 
within the corresponding radius. As expected from the ra- 
dial profiles, the gas fraction is an increasing function of 
scale radius (i.e. decreasing overdensity) . Note that, for a 
fixed enclosed mass, the variation in enclosed gas fraction is 
relatively modest: e.g. at Ma ~ lO^'* h~^ M0 it varies from 
0.07 for A = 2500 to 0.10 for A = 200. This is much less 
than the variation seen if a fixed overdensity is used to mea- 
sure the mass, e.g. for M500 = 10^* /i~^ M©, the enclosed 
gas fraction increases from 0.04 to 0.11 as the overdensity 
drops from 2500 to 200. 

In Figures [3l and [4l we fit the hot gas fraction scaling 




Figure 3. The cumulative hot gas fraction within rsoo as a func- 
tion of total mass. The solid lines show the best-fitting mean 
relations, as described in the text and Table [T] Extending the fit 
to the whole range of the FO data, as shown by the dot-dashed 
line, makes almost no difference to the fit. The dashed line shows 
the mean baryon fraction in the simulation. 
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Figure 4. The baryon fraction within three difi'erent radii for 
the FO run: r200 (upper), r^oo (middle) and r2500 (lower). For 
clarity, the upper and lower data points have been shifted by 0.5 
dex. The solid lines show the best-fitting mean relations and the 
dashed lines the best-fitting mean relation for the rsoo data. 



relations with models of the form 

logio / = logio /o + s(m - log (1 + exp(C^))/C), (3) 

where /i = logj^g(M/Mo) and /o, Mg, (" and s are fitting 
parameters. In log space, this represents a line of constant 
slope, s, at masses well below Mo, bending over to a constant 
value of /o at high masses. For low enclosed overdensities we 
would expect /o to tend towards the universal baryon frac- 
tion of 0.18, although we do not impose this as a constraint. 
^ is a parameter that controls the abruptness of the transi- 
tion between the two regimes. The data are not always suf- 
ficient to independently constrain all the parameters, and in 
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Table 1. Model parameter fits to the hot gas fractions as a func- 
tion of mass, as described by Equation |3] with a fixed value of 
f = 4. Here m = logjQ(Mo). Typical 1-sigma errors in /o, m and 
s are 0.005, 0.2 and 0.04, respectively. For the GO run there is no 
discernible mass-dependence in the hot gas fractions and so only 
the mean value is recorded in column /q. The final column gives 
the root-mean-square scatter in dex of the data points about the 
best-fit line. 



Model 


Overdensity 


/o 


m 


s 


a 


GO 


2500 


0.158 






0.036 




500 


0.161 






0.018 




200 


0.163 






0.013 




500\2500 


0.164 






0.026 




200\500 


0.166 






0.033 


PC 


2500 


0.103 


14.15 


0.590 


0.048 




500 


0.134 


14.26 


0.519 


0.041 




200 


0.150 


14.43 


0.363 


0.027 




500\2500 


0.161 


14.28 


0.393 


0.037 




200\500 


0.180 


14.10 


0.263 


0.044 


FO 


2500 


0.126 


14.25 


0.650 


0.124 




500 


0.143 


14.26 


0.552 


0.061 




200 


0.148 


14.21 


0.472 


0.058 




500\2500 


0.163 


14.13 


0.492 


0.076 




200\500 


0.173 


13.72 


0.512 


0.066 



particular C,: for that reason we use a fixed value of C, when 
recording our fits. The best-fit models are shown as solid 
lines in the figures, and the parameters are listed in Table [T] 
along with the scatter about the best-fit relation. 

Figure [S] shows the cumulative hot gas fractions in the 
FO run within r25oo and rsoo as a function of mass. The 
shaded region shows the 1-sigma spread about the mean re- 
lation from the simulation and the points are observational 
data from the sources listed in the figure caption. Concen- 
trating first on the upper panel, it is apparent that there are 
systematic differences in the reported hot-gas fractions that 
can not be attributed s olely to statistical error. In particular, 
the lAllen et al.l (|2008l ) data lie significantly above those of 
the other studies. This may be because they concentrated 
on regular, CC clusters that are likely to lie up the upper 
edge of the distribution. Even given these observational in- 
consistencies, however, the simulations show a much greater 
variation in gas fraction with M2500 than do the observa- 
tions. It would seem that in the FO run (and the PC run is 
similar) we have ejected too much material from within this 
radius in small clusters, and too little in large ones. 

Moving out to rsoo, as shown in the lower panel of the 
figure, the simulations and observations are in better agree- 
ment. There is again a suggestion that the observational 
data would prefer higher gas fractions than the simulations 
below a cluster mass of of 5 x lO^^/i'^M©, but the scat- 
ter in the observational measurements is large. At higher 
masses, the two show a similar trend of increasing gas frac- 
tion with cluster mass, although the simulated values are 
perhaps slightly too high. This is presumably because the 
mean baryon fraction that we have used in the simulations, 
0.18, is higher than the current WMAP best fit value of 
0.168. 
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Figure 5. The cumulative hot gas fraction versus mass rela- 
tions for the FO run as compared to observations. The upper 
and lower panels refer to overdensities of 2500 and 500, respec- 
tively. The shaded regions are the 1-sigma spread in the sim- 
ulated clusters. T he sy m bols represent o b servat io nal data from 
San derson et al.l ll20q3[). IVikhlinin et al.| ||2006| '). lArnaud et"ar 



(20071). lAllen et al.l ( 12008'). 'Sun et al.' ('2009'') and lPratt et al. 



(201Q!). For the lSande rson et al. (2003) clusters we show one fake 
data point on the right-hand edge of the plots with typical 1- 
sigma (statistical) error bars. For the other s amples, we include 
the error bars on the pl otted points, apart from lAllen et al.l 1120081 ) 
and iPratt et al.l ll2010|) since we do not know the uncertainties in 
their measurements. 



3.2.2 Differential hot gas fractions 

The profiles of Figure[T]suggest that the differential gas frac- 
tion between radii of r25oo and rsoo niay provide a measure 
that is more independent of mass than the cumulative gas 
fractions of the previous section. 

In Figure [6] we show differential gas fractions, r25oo- 
rsoo and r5oo-''20o, for the FO model (once again, the PC 
run gives similar results). Both are higher than the equiv- 
alent cumulative measures, although the universal gas frac- 
tion is reached only for the most massive clusters {h'hm > 
3 X lO'^* /i~^ M©) at radu r > rsoo- Clearly the differen- 
tial gas fraction at larger radii, r5oo-'"20o, is more nearly 
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Figure 6. Differential hot gas fractions for the FO runs: between 
^2500 and rsoo (lower points), and between rsoo and r200 (upper 
points). The upper points have been shifted up by 0.5 dex, for 
clarity. In each case, the mass has been taken to be that at the 
outer edge of the differential range. The solid lines show the best- 
fit mean relation and the dashed lines show the best-fit relation for 
the equivalent cumulative gas fraction measure. The dash-dotted 
lines show the mean baryon fraction in the simulation. 
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Figure 7. The differential hot gas fraction in the annulus con- 
tained between r2500 and rsoo for the FO run as compared to 
observations. The shaded region is the 1-sig ma spread in the 
simula t ed clusters. The points are t aken frorn ISanderson et al.l 
||2003|). IVikhlinin et af] 120061) and ISun et al.l ||2009| ). For the 
ISanderson et al.l 1 120031) clusters we show one fake data point on 
the right-hand edge of the plot with typical 1-sigma (statistical) 
error bars. For the other samples, we draw the error bars on the 
plotted points. 



constant and so provides the more accurate probe of cos- 
mology, but observationally the inner annulus, t-25oo~'"500 5 
provides a compromise between eliminating the depleted in- 
ner region and having enough counts to enable a reliable 
X-ray dete rmination of the gas de n sity. The observationa l 
data from ISanderson et al.l (120031 ), IVikhlinin erall (120061 ) 



Table 2. As for Table [T] but for the baryon fractions rather than 
the hot gas fractions. (^ = is fixed at 4 for the PC and 8 for the 
FO run. For the GO run there is no star formation and so the 
values are the same as reported in Table [T] 

Model Overdensity /o m s cr 

PC 



FO 



2500 


0.141 


14.32 


0.250 


0.061 


500 


0.161 


14.44 


0.307 


0.037 


200 


0.168 


14.49 


0.271 


0.030 


500\2500 


0.174 


14.29 


0.362 


0.035 


200\500 


0.188 


14.10 


0.269 


0.042 


2500 


0.156 


14.92 


0.251 


0.069 


500 


0.150 


14.48 


0.308 


0.052 


200 


0.157 


14.41 


0.291 


0.044 


500\2500 


0.167 


14.25 


0.357 


0.053 


200\500 


0.181 


13.94 


0.333 


0.058 



and lSun et al. | l|2009l). plotted in Fig ure[7l present a confused 



and loun et al. | l | zUUijll. plotted m f ig ure|Vj present a contused 
picture. The lSanderson et al.l (|2003l ) data broadly mimic the 
simulations, but the statistical scatter in their data is very 
large. Vikhlinin ot al. (2006) report the smallest error bars 
for their data and find differential gas fractions that increase 
strongly with mass, but which lie b elow the si r nulate d val- 
ues; whereas, at lower masses, the ISun et al.l (|2009J) data 
seem to require hot gas fractions that are decreasing, or 
at best flat, as a function of mass. We conclude that these 
differential measurements are not yet sufficiently robust to 
provide useful constraints. 

It is interesting to note that both the PC and FO models 
have higher differential gas fractions at large radii than does 
the GO model. The injection of entropy has removed gas 
from the cores of the clusters and pushed it out to larger 
radii, between rsoo and r2oo. In a steady-state, the higher 
entropy in these runs would ensure that they have a lower 
gas density than in the GO model: we conclude that on large 
scales, although still within the virial radius, the clusters are 
not in dynamical equilibrium. 



3. 2. 3 Baryon fractions 

The baryon fractions are also well-fit by the model given in 
Equation |3] with parameters as listed in Table [2] Because 
clusters are large systems with deep potential wells, it is 
often stated that they should enclose a representative sample 
of the Universe, and in particular that they should contain 
the Universal fraction of baryons. Indeed, we find this is 
approximately true, with only a small baryon deficit within 

r200. 

Observationally, the baryon fraction is hard to deter- 
mine because of the difficulty in measuring the contribu- 
tion from dwarf galaxies and from intracluster light (stars 
that have been stripped from galaxies). This latter com- 
ponent may compris e as much as 40 percent of the total 
light of the cluster ([Bernstein et al.l Il995l: iGonzalez et al.l 



2000l: iFeldmeier et al.l |2002|. J2004iriGonzalez et al.l l2005l : 
Krick et al.ll2006l : IZibetti et albooa ). Except for the bright- 
est X-ray clusters, the measurement of total mass is also 

pro blematic. 

iGonzalez et al.l ([20071, hereafter GZZ07), in a sample 
of 12 groups and clusters spanning a wide mass range, find 
that the baryon fraction is independent of mass and averages 
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Figure 8. The cumulative baryon fraction versus mass. The solid 
and dotted lines show the mean relations (a ppro ximately inde- 
pende nt of mass) from [Gonzalez et al.l 1120071 ) and iLagana et al.l 
1120081). respectively ; the dash-dotted line shows that from 
iGiodini et al.l 1120091 ). The red dashed line is the universal mean. 



to 0.133 within r.'inn. lLagana et al.l ([2003, hereafter LLA08) 
with a smaller sample of 5 high-mass clusters find a slightly 
lower value of 0.123 (the mean of their quoted numbers, 
weighte d by the inv erse square of their errors). On the other 
hand iGiodini et al.. (2009, hereafter GPFO^, in a sa mple of 
41 clu sters drawn fromlVikhlinin et al.l (|2006l ). lArnaud et al.l 
(J2007l ) and ISun et al.l (J2009l ). find that the baryon fraction 
is a slowly-increasing function of mass. The baryon fractions 
for our FO and PC clusters are compared to these observa- 
tions in Figure [S] 

The observations and the simulations approximately 
agree for cluster masses of 1-3 x lO^^^/i"^ Mq. However the 
simulations show a strong variation with cluster mass, even 
more so than that of GPF09. This difference is attributable 
mainly to differences in star formation, as described below. 
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Figure 9. The cumulative stellar fraction versus mass. The upper 
panel shows the stellar fractions within rsoo; the black dots are 
observed clusters from LLA08, the solid line shows the observed 
relation from GZZ07, and the dotted line that from GPF09. The 
lower panel is the stellar fraction w ithin r200 , with the solid line 
showing the observed relation from I AndreonI 1120101 ). 



3.2.4 Stellar fractions 



Both our PC and our FO models have much more mod- 
est st ar formation than do many previous simulations 



2OIOI : 



iBorgani et aP '2004'; 'Et tori et all I2OO6I; iKav et al 



2007al; iNagai et al. 2007; Dav e et all 120081 ; iFabian et al 



Puchwein et al.ll2010l ). In particular, the FO run takes 



its star-formation rate from the highly-successful L-Galaxies 
semi-analytic model. The mean stellar fraction in our high 
mass clusters in the FO run agrees well with the observa- 
tions of both GZZ07 and LLA08; the results of GPF09 are 
slightly higher. However, we do not find such a strong in- 
crease in stellar fraction in low-mass clusters as is seen in 
both GZZ07 and GPF09. The upper panel in Figure |9] shows 
the stellar mass-fraction within 7-500 as a function of mass for 
both the PC and the FO runs, with the trend from GZZ07 
shown as a solid line and that from GPF09 as a dotted 
line. This comparison suggests that we considerably under- 
estimate star formation in groups. T he lower panel s hows 
a similar plot for r2oo with data from lAndreonI l|2010l ). He 



finds lower stellar fractions but a similar steep increase with 
decreasing mass. 

Before dismissing our FO model as unrealistic, however, 
we note the following; 



• As pointed out bv lBalogh et al.1 l|2008l ). the GZZ07 data 
is incompatible with any model that forms galaxies via hier- 
archical mergers unless there is an unreasonably large star- 
formation rate in groups at late times. 

• The L-Galaxies model produces correlation functions 
for the galaxy distribution that are consistent with obser- 
vation s with no evidence for a su ppression at small separa- 
tions (JKitzbichler fc WhitdboOSl ). It is difficult to reconcile 
this with the need to greatly increase the stellar fraction in 
groups. 

The other simulations mentioned above also predict a 
slow variation of stellar fraction with mass, although they 
mostly have stellar fractions that are higher than ours, 
agreeing with the observations on group scales but having 



Baryon fractions in clusters 9 



stellar fractions that are too hish on clu ster scales. For ex- 
ample, the clusters of lEttori et alj (|2006l ) have a stellar frac- 
tion of about 0.05 within rsoo. The equivalent fraction for 
massive clusters {kT^i > 5keV) in our own runs is 0.02 for 
PC and 0.013 for FO. Q 

We conclude that, although the observations are not 
yet sufficiently robust, stellar mass fractions provide an im- 
portant test of, and discriminant between, different galaxy 
formation models. 



3.2.5 Scatter in the scaling relations 

In this section we investigate why some clusters have slightly 
more hot gas, and others slightly less, than other clusters 
of the same mass. Our purpose in doing this is two-fold: 
firstly to understand the physical reason for this scatter, 
and secondly to suggest corrections that can be applied to 
the observations to better allow gas fraction to be used as a 
probe of cosmology. 

First note that, as is evident from Figure O the mean 
gas fractions in the GO run are independent of mass. How- 
ever, there is scatter about the mean gas fraction that we 
might hope to relate somehow to the physical properties of 
the cluster. 

We have checked for correlations of the scatter with ev- 
ery conceivable physical quantity (including substructure, 
merging history, angular momentum, etc.) and find many 
weak correlations, but no strong one. Figure [10] shows a pos- 
itive correlation with concentration, i.e. more concentrated 
clusters have a greater gas fraction than the average within 
^■2500. Likewise, clusters that form earlier have a greater gas 
fraction than those that form later. The correlation coeffi- 
cients for these two relations are 0.29 and —0.23, respec- 
tively. These may be two aspects of the same relation as 
concentration shows a negative correlation with formation 
time. The appendix describes how we measure each of these 
for the clusters in our simulation. Both are correlated with 
cluster mass, but it turns out that they are more strongly 
correlated with each other. 

The physical mechanism that may drive the correlations 
seen in Figure [10] is unclear: it primarily affects the core of 
the cluster as the correlations get weaker if one measures the 
gas fraction within larger radii. It may be that the degree 
of gravitational prc-hcating is greater in systems that form 
later (see, e.g.. lMo et al.ll2005h . 

Other quantities that we have tested include the halo 
angular momentum, merger history and substructure. Once 
the primary correlation of gas fraction with mass is removed, 
none of these show any correlation with the residual gas 
fraction. 

Observationally, of course, an excess core gas fraction 
is associated with an increase in X-ray luminosity. Thus, as 
shown in Figure [TT] the excess luminosity of a cluster above 
the mean Lx-Tsi relation is correlated with the presence of 
excess gas in the cluster. The correlation coefficients in this 



^ lEttori et al.l 1120061 ) choose to give stellar fractions in terms of 
the mean baryon fraction, Y = Mstar/Aftotal//b- The stellar frac- 
tions quoted in their paper are thus a factor of 5.6 larger than 
those listed here. 
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Figure 10. The strongest correlations of the scatter in the mean 
gas fraction within r2500 for the GO run. The upper points show 
the cluster concentration, and the lower points the expansion fac- 
tor of the Universe at the time that the cluster had accumulated 
half its final mass. 
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Figure 11. The deviation from the mean gas fraction- mass re- 
lation plotted against the deviation from the mean luminosity- 
temperature relation. Properties are measured within a radius of 
rsoo- Shown are the ratio of the measured quantities compared 
to that of the best-fitting mean relation. 



case are 0.51 and 0.73 for the PC and for the FO run, re- 
spectively. A similar result was found for the PC run by 
IStanek et al.l (|20ld ). This correlation may serve to reduce 
the scatter from the major outliers in the gas fraction-mass 
relation. 



3.3 Evolution 

3.3.1 Profiles 

Figure [T5] shows the evolution of the cumulative gas fraction 

profiles of the 10 most massive clusters, for each of the runs. 

Looking first at the GO run in the upper plot, it can 
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Figure 13. Gas fractions within rsoo at 2: = 1 as a function 
of total mass. The sohd Unes show the best-fit mean relations. 
The black circles are hish-r edshift clusters (0.8 ^ z ^ 1.3) from 
the observational dataset of iMaughan et al.l 1 120081) . The dashed 
line illustrates the effect of imposing a typical survey flux limit of 
6.5x 10~^* ergs~^ cm~^ at 2; = 1; only clusters to the right of this 
lino would actually be observed by such a survey (the observed 
cluster slightly to the left of this line is at a lower redshift 2 Ri 0.8). 



profiles remain constant over time. Althougii our simulation 
takes its level of feedback from a semi-analytic model, nev- 
ertheless it seems to have achieved a homologous evolution. 
Thus, although the PC and FO clusters have indistin- 
guishable gas profiles at the current day, they look very dif- 
ferent in the past. This casts doubt on the use of the mea- 
sured gas fraction as a cosmological probe, but instead opens 
the possibility that it can be used to determine the nature 
of the feedback mechanism: in particular to distinguish be- 
tween early (PC) and continual (FO) heating. 



Figure 12. The evolution of the cumulative gas fraction profiles 
of the 10 most massive clusters. In each case the dot-dashed blue, 
dashed green, dotted magenta and solid black lines correspond to 
2 =1.5, 1, 0.5 and 0, respectively. The middle line of each set is 
the mean of, and the upper and lower the 1-sigma spread in, the 
gas fraction profiles. The red dashed line is the universal mean. 



be seen that the gas fraction within rsoo remains largely 
constant over time: this is to be expected for self-similar 
evolution. As the effective resolution increases (i.e. the ratio 
of the smoothing length to 7-500 decreases) so the gas fraction 
within the cluster core can be seen to be depleted, though 
still much higher than in the other two runs. 

In the PC simulation, the gas fractions at high redshift 
are much reduced over their current-day values. This is be- 
cause a large amount of energy has been injected into the 
ICM at early times, expelling gas from the clusters. Subse- 
quently, the gas falls back into the clusters as the Universe 
evolves and tends towards (but falls far short of) self-similar 
evolution. In other words, the early entropy injection be- 
comes relatively less important in more massive systems at 
late times. 

This is in contrast to the behaviour in the FO run. Here 
we have continual injection of energy so that gas fraction 



3.3.2 Scaling relations and comparison with observations 

Figure [131 shows the hot gas fractions within rsoo at a red- 
shift z = 1 for clusters in our three simulations. For com- 
parison, we also plot high-redshift systems (0.8 ^ 2 ^ 1.3) 
from a catalogue of clus t ers o bserved wit h Chandra com- 
piled by iMaughan et al.l (|2008h . Note that iMaughan et all 
(200a ) do not themselves present /gas values. To compute 
them, we first determine the total mass, M500, from the 
supplied values of Yx (where Vx is defined as the product 
of the gas mass within rsoo and the spectroscopic-like tem- 
perature in the spherical annulus 0.15 rsoo < r $; rsoo) by 
using the Ix — Msoo relation derived from the sample of 



Vikhlinin et al. 



Maughan et al. 



This is the procedure adopted by 
see their Equation 4). The gas frac- 
tion then follows upon taking the ratio of the gas mass in- 
terior to rsoo (tabulated in their paper) to the total mass. 
The errors on the observational data points in Figure [13] 
are computed using the supplied statistical errors on the 
gas mass and the core-excised temperatures. We compute 
errors on /gas in this way, rather than using the errors on 
Vx, because the gas mass and temperature are independent 
measurements. 

The different evolutionary behaviour of the gas fraction 
profiles is refiected in that of the mean gas fractions: the PC 
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clusters have significantly smaller gas fractions at early times 
than those in the FO run. On the whole, the predictions of 
the FO model provide a closer match to the observational 
data than those of the PC model, but it seems as if some 
of the observed data points lie below the lower edge of the 
FO /gas — Afsoo relation. This could be because the mean 
cosmic baryon fraction in our simulations is higher than that 
measured by the WMAP satellite. If we were to repeat our 
simulations with the measured value of /b, we would expect 
all relations in Figure[T3]to be shifted downwards, improving 
the agreement between our FO model and the observations. 

However, we note that the observational mass estimates 
may be lower than the true mass, because they are de- 
rived from a Ix — M500 relation that was calibrated using 
clusters with hydrostatic mass estimates. Simulations have 
shown that the assumption of hydrostatic equilibrium can 
bias such mass estim ates low by ~ 10 — 20% (Rasia ot al.l 
2006; Kav ct al. 2007b; Na gai et alJl2007l:lBurns et al.[l2008l : 
Piffaretti fc Valdarnini. ,20081 : iMeneghetti et al.l 120091 ). be- 
cause of additional pressure support provided by subsonic 
bulk motions in the ICM and/or non-thermal components. 
This would imply a small systematic over-estimate of the gas 
fraction, so the observational data points in Figure fTSl should 
be shifted downwards and to the right. Another potential 
source of systematic erro r is that the mas s es of high-redshift 
clusters in the sample of iMauehan et al.l (|200a) were deter- 
mined by assuming self-similar evolution of the Ix — M500 
relation. 

It is also important to consider the effect of source selec- 
tion on our results. Observational cluster selection is based 
on X-ray flux, so may be biased towards systems with higher 
baryon fractions, particularly at high redshift. It is not pos- 
sible t o quantify this e f fect p recisely using the archival sam- 
ple of iMaughan et alJ ()2008l ') since their selection function 
is unknown. A simple way of estimating the impact of se- 
lection effects on our findings is to ask the question: given 
a typical flux-limited survey, which of our simulated clus- 
ters would actually be observed? For illustrative purposes, 
we choose a flux limit of 6. 5 x 10~^^ e rgs~^ cm~'^, equal to 
that of the WARPS survey (|Horner et al. 20 8) fro m which 
many of the objects in the iMaughan et al.l (|2008l ) sample 
are drawn. The effect of imposing this flux limit at z = 1 
is shown by the dashed line in Figure 1131 only objects to 
the right of this line would be observed by a WARPS-like 
survey . Note that one of the clusters from IMaughan et al.l 
(|2008l ) lies slightly to the left of this line; this is because it is 
at a lower redshift, z ~ 0.8. The important point to note is 
that, for both the PC and FO runs, there is only a narrow 
mass range where the bias is signiflcant, with most clusters 
in the two samples remaining unaffected. Even if we took 
a much higher flux limit, it would require greatly increased 
scatter about the PC /gas — Afsoo relation for consistency 
with the observations, which is not intrinsic to the model. 

The best-fit parameters to the scaling relation of Equa- 
tion [3] are shown as a function of redshift in Figure 1141 In 
the case of the GO run, we fit only the mean value of the 
gas fraction, /o, which is well-determined. For the other two 
runs, the shaded regions show the 1-sigma allowed parame- 
ter range determined from monte-carlo markov chain fitting. 
The parameters show considerable scatter, but this scatter 
is highly correlated. So while it is formally possible for both 
the PC and FO clusters to have the same value of /o at 




Figure 14. Evolution of the fitting parameters of Equation |3] 
for the gas fraction within rsoo: GO (dotted, magenta on yellow), 
PC (dash-dotted, red on magenta) and FO (dashed, blue on cyan) 
lines. The shaded region in each case shows the 1-sigma allowed 
parameter region. The lowest panel shows the rms scatter in dex 
about the best-fit relation. The solid lines in the upper three 
panels show the best fit linear relations to the redshift evolution 
of the parameters. 



z = 1, for example, the other parameters must adjust them- 
selves so as to maintain the difference in gas fraction seen 
in Figure [131 

The solid lines in the upper three panels of Figure [14] 
show straight-line fits to the evolution of each of the param- 
eters with redshift - note that we fit to log^Q /o rather than 
/o as it is the former that appears in Equation [3] These fits 
are listed in Table [3] and are used in the following section 
when comparing our model predictions with observations. 

From the straight-line fits shown in Figure [T4l the mean 
gas fraction within rsoo can be predicted for clusters of 
given mass and redshift using Equation (O. This prediction 
is compared to obser ved gas fractions from the sample of 
iMaughan et all (|2008l ) in Figure [TSj What is plotted here is 
the ratio of the observed gas fraction to the predicted one, 
so that perfect agreement would correspond to a value of 
unity, independent of redshift, but with some scatter due to 
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Table 3. Best fit straight fines to the evolution with redshift of 
the gas fraction model parameters shown in Figure [141 These fits 
take the form p = po + *p2;, where p is the parameter, po its value 
at 2 = 0, and Sp the slope of the relation with redshift. 



Parameter Model 



PO 



logio /o 



14.23 


0.00 


14.30 


-0.44 


0.53 


0.08 


0.51 


0.01 


-0.870 


-0.065 


-0.834 


0.022 



Table 4. Best fit straight lines to the observed versus predicted 
hot gas ratios seen in Figure [TSl in linear and in log3/-log(l -f 
z) space. The allowed 1-sigma parameter ranges are calculated 
assuming that the expected variance about the best fit is equal 
to the observed one. The scatter is the root-mean square scatter 
about the best-fit line after allowing for the observational errors. 
In the log-log plots, the scatter is expressed in dex. 



Model 



const 



slope 



scatter 



PC linear 0.960 ± 0.022 

log 0.126 ±0.012 

FO linear 0.923 ± 0.019 



0.436 ± 0.064 0.074 
0.152 ± 0.020 0.028 
0.002 ± 0.050 0.071 



log 



-0.019 ±0.012 0.026 ±0.020 



0.034 




Fig ure 15. Ratios o f the o bserved hot gas fractions within rsoo 
from lMaughan et al.l (12003) to our model predictions. The solid 
lines show the best straight-line fits in log j/-log(l + z) space. 



the cluster-to-cluster variation and measurement error. The 
data are presented in this way since the gas fraction is a 
function of both cluster mass and redshift. The error bars 
are computed using the errors on the observed gas fractions 
(see above for details of how the se were determined from 
the data of iMaughan et al.ll2008l ), accounting for the fact 
that the error on the total mass introduces an extra uncer- 
tainty when computing the theoretical prediction for the gas 
fraction. 

It is immediately apparent that observations favour the 
FO prediction over the PC one, i.e. limited evolution in gas 
fraction since z = 1. This is shown in Table |4] where we list 
the allowable parameter ranges for straight line fits to the 
data both in linear and in logy-log(l + z) space, where y is 
the ratio of observed to predicted gas fraction within rsoo • 

In making these fits, we treat the scatter about the 
mean relation as an unknown, ascattor, independent of mass 
and redshift. The data are not good enough for a more so- 



phisticated model, and that is hkely, anyway, to make little 
difference to the fit. Given observational data, yi, and errors, 
(Ti, we estimate the scatter as 
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where ytfit are the best-fit values. We iterate to convergence 
in (Tacattor, at each stage minimising the chi-squared statistic 
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Note that the scatter about the best-fit line is, in each 
case, lower than that seen in the simulations (as shown in 
the bottom panel of Figure 1 14[). Formally, therefore, neither 
of the models is a good fit. However, it seems unlikely that 
the true scatter in /gas will be below that seen in the PC sim- 
ulation. The uncertainty in the observed /gas values is hard 
to determine, particularly at high-redshift, so it is quite pos- 
sible that the size of the error bars has been over-estimated, 
leading to an under-estimate of the intrinsic scatter. 

For the PC simulation, the slope of the observed to 
simulated gas fraction ratio is incompatible with a horizon- 
tal line with high significance. The difference between the 
best-fit values ai z — 0.1 and 2 = 1 is about 6 times the 
scatter. Even if we were to account for observational bias 
in flux-limited samples towards clusters with higher baryon 
fractions, especially at high redshift, this is simply too large 
a difference to be explained by selection effects alone (recall 
our discussion of FigurefTH)). We conclude that the PC model 
can be ruled out as a viable cause of entropy generation in 
the ICM. 

The FO simulations, on the other hand, are perfectly 
consistent with a constant ratio of approximately unity. The 
slightly lower mean hot gas fraction for the observations as 
compared to the simulations can be explained by the fact 
that the latter have a higher mean baryon fraction than the 
WMAP 7-year value.5 

We note that the analysis of Ettori et al.l (120091) has 
many clusters in common with iMaughan et al.l (|2008h , but 
lists total masses and gas fractions that are often in disagree- 
ment. We are not certain why this is but note that there are 
differences in the analysis of the data, particularly in the 
cluster outskirts. We have repe ated the analysis d escribed 
in this section with the data of lEttori et al.l 1 20091 ) but the 
data are much less constraining, principally because they 
quote much larger error bars. Nevertheless, one should bear 
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Figure 16. The relationship between the slope in log-space of the 
observed/simulated gas fraction ratio as a function of redshift and 
the value of f^A- The upper, red curve is for the PC run and the 
lower, blue curve for the FO run. The shaded regions show the 
formal 1-sigma confidence regions. 



in mind that systematic errors in the assumptions made in 
the data analysis could be degenerate with difTerences in the 
simulated ICM physics. 

Looking at the problem in reverse, one could ask what 
errors could be introduced into the determination of cos- 
mological parameters by the choice of an incorrect physical 
model for the evolution of the ICM. It is not possible to 
make a precise prediction for this using the current simula- 
tions as we only have access to a single realisation with a 
particular set of cosmological parameters. Nevertheless, fix- 
ing the simulated clusters to be the same, Figure (TH] shows 
the effect of changing the observed cluster gas fractions in 
response to different values of Q.a (fixing Q.a + fim = !)• 
From this, it can be seen that using an incorrect physical 
model can have a dramatic effect, larger than that induced 
by changing cosmological parameters within any reasonable 
range. 

The analysis described in this section is necessarily very 
naive. A full treatment would require a detailed understand- 
ing of the selection function of observed clusters, modelling 
of the scatter in the scaling relations as a function of red- 
shift, and of the mean relations as a function of cosmological 
parameters. Nevertheless, none of this is likely to alter the 
basic conclusion that both observed clusters and those de- 
rived from our FO model show little evolution in hot gas 
fractions within rsoo out to z ~ 1, whereas the PC model 
predicts a strong decline. 



4 CONCLUSIONS 

In this paper, we investigate the baryon content of clusters 
of galaxies in simulations using a variety of physical models 
for the intracluster medium: 

• GO - gravitational heating only with no radiative cool- 
ing. The purpose of this model is to test which aspects of 



the simulation evolve in a self-similar way and to provide a 
comparison for the other two runs. 

• PC - universal preheating to 200keVcm^ at 2 = 4, 
plus radiative cooling and star formation. This represents 
widespread and early heating by objects that lie below the 
resolution limit of the simulation. 

• FO - feedback taken from a semi-analytic model, in- 
cluding heating from both supernovae and active galactic 
nuclei, but without radiative cooling. The motivation for 
this model is to test heating from a realistic galaxy popu- 
lation that matches both the luminosity function and the 
black hole mass function of the current-day Universe. 

The differential hot gas fraction profiles of clusters 
in the GO simulation are approximately constant at radii 
greater than 0.2 rsoo , lying at 90 per cent of the cosmic mean. 
In the other two simulations, the profiles rise steeply from 
a low value in the cluster core before bending over to an 
approximately constant value at large radii: for the most 
massive clusters, fcTsi > 5keV, this occurs well within rsoo. 

The cumulative hot gas fraction profiles of our clusters 
in both the PC and FO runs lie well bel ow those of the 
regula r, co ol-core (CC) cluste rs observed bv lVikhlinin et aU 
(|2006l ) and I Allen et~aLl l|2008l ). However, they provide a fair 
match to the non-cool-core (NCC) clu sters found in th e 
REXCESS representative cluster survey (" Pratt et al.ll2010l ). 

When we look at integrated gas fractions within fixed 
radii, the agreement with observations is mixed. The total 
gas fraction within r2soo shows a stronger variation with 
cluster mass in the simulations than is seen in the observa- 
tions. On the other hand, the agreement within rsoo is much 
better, at least on scales above 5 x 10^''/i~^Mq. There is 
a small offset but that can be explained by the fact that 
we adopt a mean baryon fraction in our simulations that is 
higher than the current best-fit value from WMAP, respec- 
tively 0.18 and 0.168. 

A more slowly- varying function of mass is provided by 
the differential gas fraction between r2soo and rsoo- Unfor- 
tunately, the scatter in the observational data is currently 
too large to allow any meaningful comparison with the sim- 
ulations. 

In agreement with previous work, our simulated clusters 
show a much smaller dependence of stellar fraction on mass 
than is seen in observations. Our stellar fractions within rsoo 
are about 0.013 for massive clusters, Msoo > IO^^/i'^Mq, 
in the FO run, similar to observed values. The PC run 
gives slightly higher values, 0.02, whereas previous simu- 
lations can have stellar fractions as high as 0.05 (i.e. as 
much as a third of all the baryons within the cluster turned 
into stars). On the other hand, for lower- mass clusters, 
Afsoo ~ 5 X lO'^^ /i~^Mq, our mean stellar fractions of 0.015 
(FO) and 0.03 (PC) are much lower than the observed value 
of about 0.05. We note that there is some theoretical dif- 
ficulty in understanding such a steep dependence of stellar 
fraction on cluster mass, and that the observational deter- 
mination of this mass fraction is difficult, especially in low- 
mass systems. While this should prove a fruitful line of in- 
vestigation in the future, it is probably too early to draw 
firm conclusions about the validity of the models. 

We have fitted the gas fractions as a function of mass to 
relations of the form given by Equation O with the results 
shown in Table [l] The scatter about these mean relations is 
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lowest for the GO run and significantly larger for the PC, 
and especially the FO runs. Unfortunately, the observational 
data are too poor to provide an accurate measure. 

For the GO run, we might expect that the scatter about 
the mean gas fraction-mass relations has a physical origin in 
the formation history of the clusters. Indeed, there is a weak 
correlation/anti-correlation of gas fraction with concentra- 
tion/formation time (most strongly with the expansion fac- 
tor at the time that the cluster had accumulated half its 
current mass). The strongest correlation that might be used 
observationally to correct for scatter in the gas fractions is 
that between deviation from the mean Lx-T^i relation and 
the excess gas fraction. 

Although the gas fraction profiles are very similar for 
both the PC and FO runs at z = 0, their evolution is very 
different. In the former gas is heated and expelled from the 
clusters at early times, so that the gas is depleted at high 
redshift and gradually falls back into the cluster over time. 
By contrast, the continual injection of energy in the FO run 
leads to evolution that is close to self-similar. 

The evolution of halo gas fractions can therefore be used 
as a strong discriminant between models. We compare our 
simulate d clusters with t he co mpilation of Chandra clus- 
ters from lMaughan et al.l (120081 ). The observational data are 
fully-consistent with the FO predictions and disagree with 
the evolution seen in the PC simulation with high signifi- 
cance. We need to be a little careful in interpreting these 
results as this is a highly biased sample that may well con- 
tain a disproportionate number of luminous clusters at high- 
redshift. However, the scatter in gas fraction is sufficiently 
small that, even if we only observed those clusters with high 
baryon fractions at high redshift, the disagreement between 
the observations and the PC prediction would still be sig- 
nificant. We conclude that the observations favour continual 
heating, as in our FO model, over significant preheating at 
high redshift. 

A corollary of the strong dependence of gas fraction 
evolution on the physics of entropy generation is that it be- 
comes very difficult to use the gas fraction as a probe of 
cosmology. The differences caused by uncertain gas physics 
currently swamp those caused by reasonable changes in cos- 
mological parameters. In the future, as both observational 
data for high-redshift clusters and models of the ICM im- 
prove, a joint analysis should be undertaken that consid- 
ers variations in both cosmological parameters and cluster 
physics. 

The main limitation to our present study is that the 
absence of cooling in our FO simulation leaves us unable 
to model CC clusters. Whilst that does not significantly af- 
fect the gas fractions when integrated out to 7-500, it would 
be clearly desirable to also reproduce the full range of gas 
fraction profiles at smaller radii. We are working on ways 
to introduce cooling into the FO scheme without leading to 
excess production of cooled gas in CC clusters. 
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Figure 17. The abscissa shows the actual value of r200 for each 
of the clusters, whilst the ordinate shows the value predicted from 
Af25oo E'lid Afsoo assuming an NFW profile. 



Figure 18. Cluster concentration versus the expansion factor at 
the time that the clusters have accumulated one fifth of their final 
mass. The solid line shows the best-fitting power-law correlation. 
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secondary one that foUows because low mass clusters tend 
to form at lower expansion factors than more massive ones. 



APPENDIX 

The central concentration of clusters can be determined by 
measuring the mass at two different over-densities, for exam- 
ple M2500 and Msoo- Together, these uniquely d etermine the 
parameters of the NFW profile (N avarro et al. 1997) with- 
out any need to fit the profile as a function of radiuslj To test 
the effectiveness of this procedure, we show in Figure [TTl the 
value of r2oo predicted by the NFW fit (i.e. the NFW scale 
radius, a, times the concentration, 2:200) versus the actual 
value. As can be seen, the two agree very well, confirming 
that the clusters are well-fit by the NFW profile out to this 
radius. 

Cluster concentrations are often thought to have a de- 
pendence upon cluster mass, with more massive clusters hav- 
ing lower concentrations. That is indeed the case, but we find 
a much stronger correlation with cluster formation time, as 
illustrated in Figure [TS] Here the formation time is taken to 
be the value of the expansion factor when the total mass of 
all the subhalos that will go on to make up the cluster equal 
one fifth of the final cluster mass, but other definitions give 
similar correlations. We plot expansion factor rather than 
age as this gives a more linear correlation. The results are 
shown here for the GO simulation; those for the PC and FO 
runs are very slightly different because of the contribution 
of the baryons to the total mass. 

On removing the best-fitting correlation (shown as a 
solid line in the figure), the residual concentration shows no 
dependence on mass. The correlation with mass is thus a 



^ We provide IDL routines to do this 
http : //astronomy . susx . ac . uk/'^petert/nf w . pro . 
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